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/\ • Abstract A lattice version of quantum nonlinear Schrodinger (NLS) equation is considered, 

c^ \ which has significantly simple form and fullfils most of the criteria desirable for such lattice 

variants of field models. Unlike most of the known lattice NLS, the present model belongs to 
a class which does not exhibit the usual symmetry properties. However this lack of symmetry 
itself seems to be responsible for the remarkable simplification of the relevant objects in the 
theory, such as the Lax operator, the Hamiltonian and other commuting conserved quantities 
as well as their spectrum. The model allows exact quantum solution through algebraic Bethe 
ansatz and also a straightforward and natural generalisation to the vector case, giving thus 
a new exact lattice version of the vector NLS model. A deformation representing a new 
quantum integrable system involving Tamm-Dancoff like g-boson operators is constructed. 



1. Introduction 

Among integrable systems discrete models represent a special class, interest to which 
has been revived in recent years [1-3] . In the context of quantum integrable systems, apart 
from their own right as solvable quantum lattice models, they also play an important role 
by providing lattice regularised versions of the corresponding continuum models. Thus the 
lattice nonlinear Schrodinger (NLS) [4-9], lattice sine-Gordon [7] models etc. are useful for 
finding out the exact quantum solutions of the related field models through the quantum 
inverse scattering method (QISM) [10]. Moreover, these lattice versions are often able to 
unveil the hidden algebraic structures of the original field models [1,11]. 

Ideally a candidate for such quantum integrable discrete models represented by the Lax 
operator L{n \ A) should fulfill the following basic criteria: 

i) It should satisfy exactly the quantum Yang-Baxter equation (QYBE) 

R{X - fi) L{X) ® L{p) = L{p) ® L(A) R{X - fi), (1.1) 

a) It should have the same quantum R-matrix as the corresponding field model. 

Hi) The discrete Lax operator should yield the continuum one ai A ^ : 

L{n I A) -^ 1 + AC{x, X), A being the lattice constant. 

iv) The Hamiltonian of the discrete model should he 'local' and give hack the local 

field model at the continuum limit. 
Moreover, as a desirable physical requirement the Lax operator as well as the Hamiltonian 
and the energy spectrum should he as simple as possihle. 

A discrete version of NLS model was first suggested by Ablowitz and Ladik [5] . However 
the corresponding i?-matrix , a key object in QISM, is expressed through trigonometric 
functions [12] and does not coincide with the well known rational i?-matrix of the NLS field 
model, which reads: 



R{X) = AP - inl 



/ X — in 

~r ^ I . (1.2) 

\ X — IK 



That is, criterion ii) laid down above is not satisfied. Subsequently, a different discrete NLS 
was proposed [6], constructed through the Holstein-Primakoff transformation (HPT) applied 
to an infinite dimensional irreducible representation of su{2) with the classical Lax operator 
given by 

/ c -fAA+fA^^(n)*^(n ) -z^(n)*A^K(c+ f A^^(n)*^(n)) A 
^""'^ \iA^K{c+ ^A:'tlj{nYtl){n))tl){n) c + f AA + f A27/;(n)*V(n) ^ ' ^ "^^ 



related to spin parameter s = —-^ with c = 1 and canonical Poisson brackets 

{il){n),'4){mY} = 5nni ■ 

This model is free from the earlier drawback, namely it satisfies ii), though it now fails 
to fulfill the locality criterion iv) at the quantum level. As a remedy another version of lattice 
NLS was introduced [7], represented by a Lax operator which explicitly depends on lattice 
points and may be expressed as a product l{n \ A) = L(2n | \)L{2n — 1 | A), where L{n \ A) 
is taken as (1.3) with c = 1 + ^(— 1)'^kA, that is: 

L(n I A) = 



1 + i(-l)"fi:A - lAA + f A2V'(n)*^(n) -i^(n)*A Jk(1 + \{-IYkA + f A2V'(n)*V'N) 



4V -y "" I 4' 

zAy^K(l + i(-l)'^KA + f A2^(n)*7/;(n))'i/'(n) 1 + i(-l)"KA + f AA + f A2^(n)*V;(n) 

(1.4) 
This model satisfies the required criteria, but looks involved to comply with the physical 
requirement of simplicity. Finally, in a further investigation [8], a relatively simpler model 
was proposed, where the Lax operator was given directly by (1.4). However the important 
criterion iii) is now not fulfilled and the simplification achieved is also not fully satisfactory, 
as it is evident from the form of the Lax operator and structure of the following Hamiltonian: 

4 v^^, X 8 - kA , ,4 K^ . sr-^ A ,+ , /. ,-N 

« = -3;^B*" + 4 + 83^) + (3^ + ^)Ea**a. (1-6) 

n n 

where the local density tn again has different expressions, depending whether it corresponds 
to even or odd sites. For odd n it takes the form 



tr^ = {a^{n + 2)a{n + l))~^{{a^n)a{n - l))-\a^(n + l)a{n))~^ 

{a^{n + l)asa{n - l)}{a^{n + 2)a{n + 1))"\ 

with 



(1.5a) 



/K kA kA^ 1 

and at even n sites different, though similar expressions hold for both tn and a(n) [8]. The 
energy spectrum of this model obtained through the Bethe ansatz is also rather complicated, 
and gets simplified only at the continuum limit. As far as we know, until now not many other 
proposals were invoked to improve this situation [9], particularly to achieve simpler forms of 
local conserved quantities at the lattice level. However, a completely different approach was 
formulated in ref. [4] through equivalence between NLS and spin models using intertwiners of 
quantum spaces. 

Our primary aim here is to consider a quantum integrable lattice model, which at the 
continuum limit yields the more general AKNS system [13] and as an allowed reduction the 



NLS field model. The system considered fulfills all the desirable requirements of a discrete 
quantum system listed above and most importantly, exhibits considerably simple expressions 
for the related conserved quantities at the lattice level. Indeed it satisfies the QYBE exactly 
with the same rational i?-matrix (1.2) allowing also solution via QISM, yields local Hamilto- 
nian and gives back both the Hamiltonian and the Lax operator of the NLS field model at the 
continuum limit. Moreover, it has an extremely simple structure, which induces an almost 
trivial form for the projector required for the construction of local Hamiltonians. Remark- 
ably, this projector turns out to be field-independent and symmetric. As a further relevant 
feature, our model also allows a natural vector generalisation at the lattice level. Finally 
it admits an integrable deformation involving Tamm-Dancoff type g-bosons. On the other 
hand, for achieving all these nice properties, one has to pay some price which is reflected in 
the non-hermitian nature of the physical observables at the lattice level. At the same time 
the associated Lax operator lacks the usual SU{2){SU{1, 1)) symmetry. 

We should stress here that such Lax operators with lesser symmetries were found also 
to be significant in generating a large class of quantum integrable models [11]. 

2. The classical model 

The model under scrutiny may be given through the Lax operator of the form 

j^(^| ^) = /^C + AK(/)(n)7/;^(n) -i0(n)(KA)i\ .^j) 

Its simplified structure compared to (1.3) is explicit, though due to no-conjugacy relation 
between between (p and t/', it is obviously not hermitian. Note that similar forms of L 
operators appear also in analysing descrete self-trapping systems [14] as well as integrable 
systems close to Toda lattice [15]. 

Recently, the bi-hamiltonian structure of of the classical system corresponding to (2.1) 
has been determined and its complete integrability has been established [16] rigorously 
through the explicit construction of action-angle variables using the r-matrix approach. Re- 
call that at the classical limit the QYBE (1.1) reduces to the classical Yang-Baxter equation 

{L{n I ®, L{m \ v) } = [ r{C, v) , H^ \ () (g) L{m \ if) ](5nm • 

For the present model, the quantum i?- matrix is given by (1.2) and is related to its classical 
counterpart r by: 

^Pi?(C)=I-mr(C), r(C) = |p. 

Now to show the transition of the Lax operator to that of the continuum model one 
should put C = 1 + ^A-^ and ^(n) — > i^/Aijj{n),(j){n) -^ — zV^0(n), which introducing 



ip{n) = ^ /^" 'ip{x) dx , and a similar expression for 0, would yield from (2.1) : L{n 



A JXr^ 

C) = l + A/:(x,A) + 0(A2). 

£(a;, A) is the Lax operator of the corresponding field model, given by 

It may be easily checked that the conserved quantities associated with this system are the 
same as those of the AKNS system [13]; moreover, since their Poisson structures coincide, 
one may conclude that the two systems are equivalent. In fact through a simple gauge 
transformation 

C^hCh-^ + Kh-^ , /i = e-^^^ (2.3), 

this ^-operator can be changed into the standard Lax operator of continuum NLS: 

i^ K2(p\ .A 

•A 



£{x,X)={ \\ l]=i-a^ + n-^(l)a+ + K-^tlja-, (2.4) 



restoring the unitary symmetry, since as is well known, the AKNS system allows the reduction 

3. Quantum model 

Recently, more general forms of L-operator of discrete quantum integrable models corre- 
sponding to standard /^-matrices have been proposed [11]. Such class of L operators associ- 
ated with the rational i?-matrix (1.2) and satisfying the QYBE may be given by the following 
expression which clearly lacks the unitary symmetry: 

K^ Ks + i^^Kj^ ^^-^^ 

where K operators satisfy the algebra 

[K+, K.] = {KiK4 - K2Ks), [Ki.Ks] = 

[KuK±]=±K±K2. [Ks.K±]=tK±K4, 

with K2, K4 as central elements. It may be seen that when Ki = —K^, K2 = K^ = 1, K^ = 
{K-Y , (3.2) reduces to the standard su{2) algebra and one can get back the known lattice 
NLS (1.3) through HPT. However the L operator (3.1) in general gives the possibility of 
generating other quantum integrable models wich do not exhibit such symmetry. Quantum 
Toda chain is one of the main examples [11]. It is interesting to observe that the quantum 
version of the NLS model (2.1) considered here, also falls into this class and can be obtained 
from (3.1) through the following realisation: 

Ki = /\^kH + 1, K2 = -Ak, Ks = 1, K^ = 0, K+ = i^^fni), K- = -iAy^cp, (3.3) 



where the operators -0, (p obey the canonical commutation relation [ ip{n) , 0(m) ] = -^Smn ■ 
This quantum model, represented by the Lax operator 

_fl-iXA + A^K(P{n)ilj{n) -zAK^(n)\ 
\ lAK^ipl^n) 1 / 

as descendant of the integrable 'ancestor' model (3.1), is naturally quantum integrable and 
satisfies the QYBE with the same R matrix (1.2) as the NLS field model. 

In exact analogy with the classical case, (3.4) allows transition to the Lax operator of 
the AKNS system and through allowed reduction to that of the continuum quantum NLS 
model. Indeed, the gauge transformation (2.3), being independent of the field operators, is 
clearly applicable to the quantum case as well. 

It is known [7] that the conserved quantities Ci may be obtained from the transfer matrix 
r(A) = tr [YIn -^(^ I '^) ) 5 through an expansion at a special point u, in the form: 



In what follows we use the method developed in [7,8]. The locality of the Hamiltonian 
and other conserved quantities can be achieved provided that at this special point u the 
operator L{X) be expressible both as a "direct" and an "inverse" one-dimensional projector 
[7,9]. This in turn implies the vanishing of its quantum determinant [6] detqL at this point, 
where 

detqL = tr{V- {L{X) (g) L{X + iK.))) 

— -^[(-^11-^22 + -^22-^11) — (-^21-^12 + -^12-^21)] 

with V- = |(1 — ^^ aa ® o"a) being the antisymmetriser and L = L{X), L = L{X + in). 

We observe that for Lax operator (3.4) one gets detqL = 1 — iXA, giving a single de- 
generacy point z/i = — ^. The resulting projector depends on the field operators and one 
cannot avoid the implementation of the involved procedure discovered and applied in [6-8] 
and elaborated in [9] . Fortunately however, under an irrelevant scaling of the Lax operator 
L ^ L = i^)L, which evidently does not affect the QYBE and thus can only give equivalent 
lattice models, the quantum determinant becomes 



" A2 \XiX + iK)J A2(1 + kO' 

where C = j- That is, another degeneracy point ^ = z/2 = naturally appears. 

6 



The rescaled operator L takes the form 

with N{k) = 1 + K/S?(j){k)'^{k). At the new degeneracy point ^ = z/2 = 0, it becomes 
remarkably simple as it turns into a field independent projector: 

^(0)=(J l^=V (3.5) 

The above procedure amounts essentially to choose the expansion point at A = oo. We 
emphasize that the existence of such an exceptional expansion point where the projector 
becomes field-independent is possible only due to the asymmetry of the present model. We 
note incidentally that an analogous property also holds at the classical level [16]. As a 
consequence, due to the almost trivial form of L{Q) (3.5), as we will see now, not only the 
required locality condition is satisfied, but also the derivation as well as the expression of the 
Hamiltonian and the other conserved quantities become extremely simple. 

For explicit calculations we use now L and expand around ^ = assuming periodic 
boundary conditions, dropping however the hat sign from all subsequent expressions. This 
gives 

r(0)=tr(j]L(A;|OI(c=o))=l, (3.5a) 

^r(0|c=o = r'(0)=tr5^(L(A^|0---i^'(fc|0---i^(l|0)|(e=o). 

k 



k 



(3.56) 



A 

k 



In a similar way one gets 



r"(O) = 2|i^5^Ar(^)A^(A;) + /.5^0(^ + l)^(fc)j, (3.5c) 

V i>k k / 

where a factor 2 appears due to the identity 

(■■•L(z)...L'(fc)...)'k=o=(---^'(^)---i(fc)---)'k=o 
valid since L"{^) = 0. Continuing further we get 

^"'(°) = l^i^ 5Z Nit)NU)N{k)+K{ Y. N{i)<Pik + l)^{;{k) + 

i>j>k i,k{i=lk,^k+l) 



J2<i>ik+imk-m, 



(3.5rf) 



etc. Notice that the conserved quantities Ck may be given through the above expressions 
(3.5) in the foUowing form 

Ci = -(/o^r(O)' I (? = 0) = -r(0)-V(0) 
C2 = ^i^ogriOr I (^ = 0) = ^[r{Or'r"{0) - {r{0)-'r\0)f] 



^3 = T^(/oM0)'"l(e = 0) 

= j^[2{r-\'{0)f + r-\"'{0) - 2{r-\'{0)){r-\"{0)) - (r-V"(0))(r- V'(0))], 

where r~^= r~^(0) Inserting now the expressions (3.5) one finally gets the required observ- 
ables 

7V = Ci = -^ViV(A;) (3.6a) 

k 

as the 'Number ' operator, 

k k 

as the 'Momentum' operator and 

H = C^ = Y,hk = ^ J2i<P{k + imk - 1) - {N{k) + N{k + l))4>{k + imk) + 

k k 

{3KAY^N{kf). 



(3.6c) 



as the Hamiltonian of the system. 

It may be noted that the above conserved quantities are not symmetric in and t(j, which 
is a consequence of the asymmetry of the Lax operator. On the other hand their locality is 
explicit and it is interesting to observe that even though expressions (3.5) given through 
expansion of r(^) were all nonlocal, in the corresponding conserved quantities (3.6) all such 
nonlocal terms get cancelled among themselves leaving only the local ones, as it occurs also 
in the classical case. We stress again that the evident simplicity of expressions (3.6 a,b,c) for 
the conserved quantities, is the most prominent feature of the present model. 

The transition of these conserved quantities to those of the NLS field model is easily 
achieved at the continuum limit by taking 

^=^^ E(^(^) - 1)) l(A-o)= I dx<j){x)ij{x), (3.7a) 



P = 2{J2(.Pk + 1^)) liA^o,^ J dx{4,,i- - H:), (3.76) 



^ (3.7c) 

c?a;(0^?/^^ + K(</)V^)^), 



3kA2-I(^^°) 



with the standard assumption of vanishing boundary condition. It is worth remarking that 
the continuous conserved quantities (3.7) of the AKNS type system are now symmetric in 
and ip, which aUows therefore the reduction = ^^ yielding the known expressions for the 
NLS field model. 

The evident closeness between the conserved quantities (3.6) of the lattice version with 
those of (3.7) related to the continuum model is a noticable feature of the present model. 
For solving the eigenvalue problem for the Hamiltonian of the discrete model exactly, we 
go along the well established steps [8] of algebraic Bethe ansatz forming the basic tool of 

QISM [10]. Defining the monodromy matrix as T(A) = I ^/\\ n/'Xw ~ ^k -^(^ I ^) ^^ 
get the expression for transfer matrix as r(A) = tr{T{X)) = A{X) + -D(A), which generates 
the conserved quantities, while i?(A),C(A) acts as 'creation' and ' annihilation' operators, 
respectively. The n-particle eigenstates may be defined as | n >— Yli ^(^i) I ^ ^ with the 
property of the 'vacuum' : 

C(A) I >= 0, ^(A) I >= a(A)^ | >, D{X) \ >= d{X)^ | > . 

The QYBE for the monodromy matrix is given again by equation (1.1), with L operators 
replaced by the corresponding T operators. In elementwise form this equation yields the 
'commutation' relations 

[^(A),^(,,)l = lD{X),D{^L)] = lB{\),B(p)] = [C(A),C(A)] = 

^(^'^C^' - ^^f'-^*^' - Sj^^'^'-^f' (3.8) 



with 



H>'^l^) = -. -, c(A,^) 



X — n — iK X — n — in 

The eigenvalues of r(A) giving the physical observables may be obtained by using the com- 
mutation relations (3.8) between A, B and D, B and the properties of the 'vacuum' stated 
above. Skipping out the details we present only the main results as follows. 

r(A) |n>=^(A,{Aj}) |n>, (3.9) 

9 



where 

f(A,{A,})^a(A)-n^(^+''(A)"n^(V^ (3.10) 

Note that the form of eigenvalues (3.10) is obtained provided the parameters A, satisfy the 
condition [10] 

(^)"=n7r4T' ^'^■=i'---'- (3-11) 

For the present case one obtains a(f) = (1 + ^), d{C,) = ^. We define the Hamiltonians 
by 

where a~^ is included to remove irrelevant constant terms and to avoid linear combinations 
of conserved quantities. Thus we get from (3.10) Ck = ^^1^ logf(^) 1(^=0); where 

r(o=n '+;'^^rr'^ . (3,12) 



which finally yields 



Af=-Ci= n, 



n 
AC 



k=i 
2 ' 



2 ■ ^ (3.13) 



k=i 

We observe that the energy is proportional to A^, the momentum is proportional to \k and 
the number of particles is equal to the quasiparticle excitation number, as required. Note 
again that this result conserning the descrete model under consideration is much similar to 
that of the NLS field model [10] including the combinations of different conserved quantities 
to determine the momentum and the energy spectrum. However, contrary to the continuum 
case, here the values of A^'s are not arbitrary and should be determined from the equations 
(3.11), which for the present model reads: 



(l-.A,A)^ = n^^-4^, ^,J = l,---,n (3.14) 

1^3 



Aj — Afc + m 



We should emphasise that the energy spectrum of this model obtained above and the related 
constraints on Afc are indeed extremely simple. 

4. Vector generalisation of the model 

10 



It is interesting to observe that the models violating the SU{2) type symmetry, proposed 
in [11] can be easily generalised for the gl{N) case. Out of such generalised system one might 
then construct quantum integrable models, like multi-component Toda chain, vector NLS 
model etc., as realisations through a set of independent bosonic operators. Such generalised 
systems are given by the Lax operator 

L = Y.^K^ + T^Den + E^^^-^^-^' (4.1) 

rv 

where {eij)ki = SikSji are the generators of gl{N). It can be shown that the above L operator, 
associated with the rational (A^^ x A^^) i?-matrix : R{X) = l + i^II, where 11 = ^^^ eki ®eik^ 
satisfies the QYBE if the generators K yield the following algebra: 

[Krak, Kki] = K-Kmi, [Kki, Kik] = K+ Kf - K-K+, 

[K+, Kki] = KkiK-, [K+, Kik] = -KikK-, (4.2) 

where i^^ commute with all other generators: [Kj^ ,Kij] = and thus are central elements, 
while Ki form an abelian subalgebra. We may notice again that in general this is not a 
su{N) algebra, which however is recovered at some particular symmetric reduction. 

Different realisations of this algebra would generate through (4.1) different quantum 
integrable models, which would share the same rational i?-matrix but generically would not 
exhibit unitary symmetry. Consider now a realisation of (4.2) through a set of independent 
operators with the commutation relations [ipi, (pk] = Sik, [i^i, ipk] = [<^i, 0fc] = 0, in the form 

(4.3) 
Kij = ^,, Kji = 0,, K,j =0, 1< (ij) < N 

The corresponding Lax operator (4.1) will then read 

This Lax operator, which yields a quantum integrable lattice model, gives the vector NLS 
model [17,18] at the continuum limit and is a natural generalisation of (3.4) to the vector 
case. The associated i?-matrix also coincides with that of the field model [18]. Thus (4.4) 
is related to the Lax operator of a new exactly integrable lattice version of the vector NLS 
model. The corresponding classical system has been considered in [16]. 

5. A novel quantum integrable Tamm-Dancoff g-bosonic model 

11 



A number of lattice models involving g-oscillators , which are integrable at the quantum 
level have already been discovered [11,19]. Most of these models are related to the quantum 
group structures associated with the trigonometric R matrix, which forms a separate class 
entirely different from the NLS model with rational R matrix (1.2). We present here an 
integrable deformation of the discrete NLS model (3.4), which involves Tamm-Dancoff (TD) 
[20] type g-boson operator, but at the same time is related to a rational i?-matix. 

It has been shown in [11] that for a 'Symmetry breaking' transformation [21]: i?^- (A) -^ 
Qid{j-k) Ji>pj(^X)^ of the original R-matrix (1.2), where 9 is some constant parameter, the 
algebra (3.1) of K operators is also deformed in an interesting way. We find a realisation of 
this deformed algebra through TD type g-bosonic operators b, c, N: 

[b,N] = b, [c,N] = -c, bc-qcb = q^, (5.1) 

where b and c are not in general hermitian conjugate of each other. Here we have introduced 
the parameter q = e*^^. One may compare the above TD type (/-deformed bosons with the 
standard (/-oscillator algebra [22]: [a,N] = a, [a^,N] = —a"!", ,aa^ — qa^a = q~^ . The 
algebraic relations (5.1) yield the Lax operator of the corresponding model as 

"" ^^^ - I -^ h f{N) ) ' ('•'^' 

where f{N) = q2(^~2\ 

Note that it represents a quantum integrable system, which satisfies the QYBE with the 
deformed rational R matrix 



R\X) 



/ X — in \ 

—in Ag~2 

Ag2 —in 
\ X — in / 



(5.3) 



Evidently at q -^ 1 one recovers from (5.3) the Lax operator (3.4) of our discrete NLS and 
also gets R'^ ^ R 3iS in (1.2). 

There is a simple mapping from such TD-deformed operators to the operators of the 
original lattice NLS model as b = f{N)'ip ,c = f{N)(f) recovering the canonical relation 
[iIj, (p] — 1; accordingly the Lax operator (5.2) is mapped into (3.4) by: 



L'J(n) = g5(^(-)-5)L(n) 



NLS- 



Hence, this Tamm-Dancoff type deformed bosonic system represents a new quantum inte- 
grable lattice model which can be solved through the algebraic Bethe ansatz using the results 
reported in sect. 3. 

12 



6. Concluding remarks 

It might be wortliwliile to summarize here the main results obtained in this paper, 
and to stress once again its underlying ideology. We have presented a quantum integrable 
lattice model, which in general corresponds to the AKNS type system. It may also be 
considered as a lattice NLS model, that shares with the continuum field model the same 
quantum R-matrix, yields local and quite simple expressions for the conserved quantities 
including the Hamiltonian and allows to determine their spectrum by the standard quantum 
inverse scattering method. Remarkably, the present lattice NLS admits a natural vector 
generalisation giving a new exact lattice version of the vector NLS model, which is much 
simpler than all other vector generalisations available in the literature. Finally, the model 
can be easily deformed, giving rise to a Tamm-Dancoff type q-boson model, exactly solvable 
at the quantum level. The price we had to pay to achieve all the previous results was the 
breaking of unitary symmetry, which is however restored at the continuum limit. We note that 
the same advantages and the same drawback characterise the classical version investigated 
in [16]. For those who think, as we do, that the price we paid was an acceptable one, the 
ideology underlying the construction of the present model might be a fruitful one. 

After the completion of this work ref. [23,24] were brought to our notice. In [23], as a 
significant contribution, a most general form of L operator for the lattice NLS was found, 
which provides the basis for classification of all L operators related to the i?-matrix (1.2). 
The physical and mathematical properties of lattice NLS along with many other models have 
been discussed in great detail in [24]. We note that the model presented here is consistent 
with the general form of the L operator of [23] , which thus gives another basis for its validity. 
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